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C$ , Abstract. We consider nonlinear elliptic problems involving a 

nonlocal operator: the square root of the Laplacian in a bounded 
domain with zero Dirichlct boundary conditions. For positive solu- 
tions to problems with power nonlinearities, we establish existence 
and regularity results, as well as a priori estimates of Gidas-Spruck 
P I ■ type. In addition, among other results, we prove a symmetry the- 

| orem of Gidas-Ni-Nirenberg type. 

a 

1. Introduction 

^ ! This paper is concerned with the study of positive solutions to non- 

| linear problems involving a nonlocal positive operator: the square root 

CN ■ of the Laplacian in a bounded domain with zero Dirichlet boundary 

conditions. We look for solutions to the nonlinear problem 

in 
o 

On 

{ u = ondtt, (1.1) 

in n, 




f(u) in tt, 



^ . where Q is a smooth bounded domain of M n and Ai/ 2 stands for the 

square root of the Laplacian operator —A in with zero Dirichlet 
boundary values on dQ. 

To define A 1 / 2 , let {A^, ipk}kLi be the eigenvalues and corresponding 
eigenfunctions of the Laplacian operator —A in Q with zero Dirichlet 
boundary values on dQ, 

-Aip k = \ k ip k in Q, 
ifk = on Q, 
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normalized by H^Hl 2 ^) = 1. The square root of the Dirichlet Lapla- 
cian, that we denote by Ai/ 2 , is given by 

oo oo 

u = ^2 Cki Pk 1 — > Ai/ 2 u = ^2 c k><l /2 Vk, (1.2) 

k=l k=l 

which clearly maps Hq(Q) = {u = J2T=i c Wk \ Y^k=\ -^ c fc < oo} into 
L-(il). ^ 

The fractions of the Laplacian, such as the previous square root Ai/ 2 , 
are the infinitesimal generators of Levy stable diffusion processes and 
appear in anomalous diffusions in plasmas, flames propagation and 
chemical reactions in liquids, population dynamics, geophysical fluid 
dynamics, and American options in finance. 

Essential to the results in this paper is to realize the nonlocal opera- 
tor A1/2 in the following local manner. Given a function u defined in Q, 
we consider its harmonic extension v in the cylinder C := Q x (0, oo), 
with v vanishing on the lateral boundary OlC := dVt x [0, oo). Then, 
A1/2 is given by the Dirichlet to Neumann map on f2, u i— > §^|nx{o} ; 
of such harmonic extension in the cylinder. In this way, we transform 
problem (11.11) to a local problem in one more dimension. By studying 
this problem with classical local techniques, we establish existence of 
positive solutions for problems with subcritical power nonlinearities, 
regularity and an L°°-estimate of Brezis-Kato type for weak solutions, 
a priori estimates of Gidas-Spruck type, and a nonlinear Liouville type 
result for the square root of the Laplacian in the half-space. We also 
obtain a symmetry theorem of Gidas-Ni-Nirenberg type. 

The analogue problem to ( 11.11) for the Laplacian has been investi- 
gated widely in the last decades. This is the problem 

—Au = f(u) in Q, 

u = on dn, (1.3) 

u > in Q ; 

see [21] and references therein. Considering the minimization problem 
min{||u||jji(n) | ||w||ip+i(n) = 1}, one obtains a positive solution of (11.31) 
in the case f(u) = u p ,l<p< since the Sobolev embedding is 
compact. Ambrosetti and Rabinowitz [1] introduced the mountain pass 
theorem to study problem (ll.3p for more general subcritical nonlineari- 

ties. Instead, for f[u) = u n ~ 2 } Pohozaev identity leads to nonexistence 
to (11. 3p if Q is star-shaped. In contrast, Brezis and Nirenberg [1] showed 
that the nonexistence of solution may be reverted by adding a small 
linear perturbation to the critical power nonlinearity. 
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For the square root Au2 of the Laplacian, we derive the following 
result on existence of positive solutions to problem (11.11) . 

Theorem 1.1. Let n > 1 be an integer and 2" = when n > 2. 
Suppose that Q is a smooth bounded domain in M. n and f(u) = u p . 
Assume that 1 < p < 2^ — 1 = ^j- if n > 2, orl<p<ooifn = l. 

Then, problem (II. ip admits at least one solution. This solution (as 
well as every weak solution) belongs to C 2 ' a (Q) for some < a < 1. 

As mentioned before, we realize problem (II. ip through a local prob- 
lem in one more dimension by a Dirichlet to Neumann map. This 
provides a variational structure to the problem, and we study its cor- 
responding minimization problem. Here the Sobolev trace embedding 
comes into play, and its critical exponent 2" = -^j, n > 2, is the power 
appearing in Theorem 11.11 We call p critical (respectively, subcritical 
or supercritical) when p = 2^ — 1 = j^j- (respectively, p < 2" — 1 or 
p > 2" — 1). In the subcritical case of Theorem 11.11 the compactness 
of the Sobolev trace embedding in bounded domains leads to the exis- 
tence of solution. Its regularity will be consequence of further results 
presented later in this introduction. 

Remark 1.2. In [26] the second author J. Tan establishes the non- 
existence of classical solutions to (II. ip with f(u) = u p in star-shaped 
domains for the critical and supercritical cases. In addition, an exis- 
tence result of Brezis-Nirenberg type jl] for f(u) = u v + fiu, // > 0, is 
also established. 

Gidas and Spruck [H] established a priori estimates for positive solu- 
tions of problem (ll.3p when f{u) = u v and p < Its proof involves 
the method of blow-up combined with two important ingredients: non- 
linear Liouville type results in all space and in a half-space. The proofs 
of such Liouville theorems are based on the Kelvin transform and the 
moving planes method or the moving spheres method. Here we estab- 
lish an analogue: the following a priori estimates of Gidas-Spruck type 
for solutions of problem (II. ip . 

Theorem 1.3. Let n > 2 and 2 tt = Assume that Q C JR" is a 

71—1 

smooth bounded domain and f(u) = u v , \ < p <2^ — 1 = 

Then, there exists a constant C(p, Q), which depends only onp andVL, 
such that every weak solution of (11.11) satisfies 

IMU°°(n) < C(p,tt). 

To prove this result, we combine the blow-up method and two useful 
ingredients: a nonlinear Liouville theorem for the square root of the 
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Laplacian in all of R n , and a similar one in the half-space WL with zero 
Dirichlet boundary values on cM™. The first one in the whole space 
was proved by Ou [22] using the moving planes method and by Y.Y. 
Li, M. Zhu and L. Zhang |18j . [IT] using the moving spheres method. 
Its statement is the following. 



n+l 
n— 1 ; 



Theorem 1.4. ([IB]. [22] . [17] ) For n > 2 and 1 < p < 2 tt - 1 

t/iere exists no weak solution of problem 

(1.4) 




u p in R n 



in R n . 




As we will see later, here (—A) 1 / 2 is the usual half-Laplacian in all 
of R n , and problem (11.41) is equivalent to problem Av = and v > 
in d v v = v p on cM™ +1 . The corresponding Liouville theorem for 

the square root of the Laplacian in = {x G R n \ x n > 0} was not 
available and we establish it in this paper for bounded solutions. 

Theorem 1.5. Let n > 2, 2« = and 1 < p < 2 s - 1 = ^±± Then, 

— ' n—l' r — n— 1 ' 

i/iere exisfo no bounded solution u of 

--v? mRl, 

ondRl, (1.5) 
in M™, 

where Ai/ 2 is the square root of the Laplacian in = {x n > 0} with 
zero Dirichlet boundary conditions on cM™ . 
In an equivalent way, let 

RlX 1 = {z = (x u x 2 , ■ ■ • , x n , y)\x n >0,y> 0}. 

If n > 2 and 1 < p < 2^ — 1 = 2±1^ ^ en there exists no bounded 

solution v G C^RlX 1 ) n C(M™+ 1 ) o/ 

fAw = ml^ 1 , 

u = on {x n = 0,y > 0}, 

g = «P on{x n >0,i/ = 0}, 

v>0 inRlf, 

where v is the unit outer normal to at {x n > 0, y = 0}. 

The proof of this result combines the Kelvin transform, the moving 
planes method, and a Hamiltonian identity for the half-Laplacian found 
by Cabre and Sola-Morales [5] . The result of Theorem 11.51 is still open 
without the assumption of boundedness of the solution. 
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Gidas, Ni, and Nirenberg [13] established symmetry properties for 
solutions to problem (ll.3p when / is Lipschitz continuous and Q has 
certain symmetries. The proof of these symmetry results uses the max- 
imum principle and the moving planes method. The moving planes 
method was introduced by Alexandroff to study a geometric problem, 
while in the framework of problem (11.31) was first used by Serrin. In the 
improved version of Berestycki and Nirenberg [3], it replaces the use of 
Hopf's lemma by a maximum principle in domains of small measure. 

Here we proceed in a similar manner and obtain the following sym- 
metry result of Gidas-Ni-Nirenberg type for (11. ip . 

Theorem 1.6. Assume that Q is a bounded smooth domain of W l 
which is convex in the x\ direction and symmetric with respect to the 
hyperplane {x\ = 0}. Let f be Lipschitz continuous and u be a C 2 (f2) 
solution of (II. ip . 

Then, u is symmetric with respect to X\, i.e., u(—Xi,x') = u{x\,x') 
for all (xi, x') E n. In addition, J^- < for X\ > 0. 

In particular, if Q = Br(0) is a ball, then u is radially symmetric, 
u = u(\x\) = u(r) for r = \x\, and it is decreasing, i.e., u r < for 
< r < R. 

We prove this symmetry result by using the moving planes method 
combined with the following maximum principle for the square root 
A1/2 of the Laplacian in domains of small measure (see Proposition 14.41 
for a more general statement in nonsmooth domains). 

Proposition 1.7. Assume that u G C 2 (Q) satisfies 



where Q is a smooth bounded domain in W 1 and c G L°°(Q). Then, 
there exists 5 > depending only on n and 1 1 c such that if 

\VL n {u< 0}| < 5 then u>0 inVt. 

The above maximum principle in "small" domains replaces the use of 
Hopf's lemma to prove symmetry results for A\/2 in Lipschitz domains. 
We point out that Chipot, Chlebfk, Fila, and Shafrir [9] studied a 
related problem: 



where /, g G C 1 (1R) and v is the unit outer normal. They proved 
existence, non-existence, and axial symmetry results for solutions of 
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( 11.71) . Following one of their proofs, we establish Hopf 's lemma for Ax/2, 
Lemma [4.31 below. Finally, let us mention that singular solutions and 
extremal solutions of similar problems to (11.71) have been considered 
by Davila, Dupaigne, and Montenegro [10J, [llj. 

As we mentioned, crucial to our results is that Ax/2 is a nonlocal 
operator in fl but which can be realized through a local problem in 
Q x (0, oo). To explain this, let us start with the square root of the 
Laplacian (or half-Laplacian) in R n . Let a be a bounded continuous 
function in all of R n . There is a unique harmonic extension v of u in 
the half-space R™ +1 = R™ x (0, oo). That is, 

Av = in R™ +1 = {(x, y) G R™ x (0, oo)}, 
u onR n = <9R™ +1 . 

Consider the operator T : u i— > —d y v(-, 0). Since d y v is still a harmonic 
function, if we apply the operator T twice, we obtain 

(T o T)u = dyyV \ y =o— — A x v \ y =o— —Au in R™. 

Thus, we see that the operator T mapping the Dirichlet data u to the 
Neumann data —d y v(-,0) is actually a square root of the Laplacian. 
Indeed it coincides with the usual half-Laplacian, see [16] . 

Here we introduce a new analogue extension problem in a cylinder 
C := fl x (0, oo) in one more dimension to realize (11.11) by a local 
problem in C. More precisely, we look for a function v with v (■, 0) = u 
in R n satisfying the following mixed boundary value problem in a half- 
cylinder: 

' Av = inC = Q x (0,oo), 

on d L C := dQ x [0,oo), 



£ = /(«) onfix{0}, 



v > in C 



where v is the unit outer normal to C at Q x {0}. If v satisfies (ll.8p . then 
the trace u on fix {0} of v is a solution of problem (11.11) . Indeed, since 
d y v is harmonic and also vanishes on the lateral boundary dfl x [0, oo), 
we see as before that the Dirichlet to Neumann map u t— > —d y v(-, 0) is 
the unique positive square root Ax/2 of the Dirichlet Laplacian in Q. 

The generators of Levy symmetric stable diffusion processes are the 
fractional powers of the Laplacian (— A) s in all of R™, < s < 1. 
Fractional Laplacians attract nowadays much interest in physics, bi- 
ology, finance, as well as in mathematical nonlinear analysis (see [2]). 
One of the few nonlinear results which is not recent is due to Sugi- 
tani [25J, who proved blow up results for solutions of heat equations 
d t u + (— A) s u = f(u) in R n , for all < s < 1. It is important to note 
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that the fundamental solution of the fractional heat equation has power 
decay (or heavy) tails, in contrast with the exponential decay in case 
of the classical heat equation. Levy processes have also been applied to 
model American options |2J. As recent nonlinear works for fractional 
diffusions, let us mention the following. Caffarelli and Silvestre [7] have 
given a new local realization of the fractional Laplacian (— A) s , for all 
< s < 1, through the Dirichlet-Neumann map of an appropriate de- 
generate elliptic operator. The regularity of the obstacle problem for 
the fractional powers of the Laplacian operator was proved by Silvestre 
[23] . The optimal regularity for such Signorini problem was improved 
in [BJ. Moreover, the operator (— A) s plays an important role in the 
study of the quasi-geostrophic equations in geophysical fluid dynamics; 
see the important recent paper [8] by Caffarelli and Vasseur. Cabre and 
Sola-Morales [5] studied layer solutions (solutions which are monotone 
with respect to one variable) of (— A) 1//2 w = f(u) in W 1 , where / is of 
balanced bistable type. 

To prove Theorem 11.11 in view of (11.81) being a local realization of 
(11 .11) . we consider the Sobolev space 

H l L {C) = {v G H\C) | v = a.e. on d L C = dVl x [0, oo) }, 

1 /2 

equipped with the norm ||t>|| = (J c | Vf | 2 dxdy) . Since problem (11.81) 
has variational structure, we consider its corresponding minimization 
problem 

J = inf|y \Vv(x,y)\ 2 dxdy\veHl L (C),J \v(x, 0)\ p+1 dx = 

We will prove that, for subcritical powers, there is a minimizer for this 
problem. Its trace on Q x {0} will provide with a weak solution of 

m- 

Thus, it is important to characterize the space Vq(Q) of all traces on 
Q x {0} of functions in H^ L (C). This is stated in the following result 
— which corresponds to Proposition 12.11 of next section. 

Proposition 1.8. Let Vq(Q) be the space of all traces on Q x {0} of 
functions in Hq L (C). Then, we have 

Vo(ft) := {u = tr n v\veHl L (C)} 

= {ueH 1 ' 2 {Q)\ [ ^Idx <+oo\ 

oo oo 

u G L 2 {Q) \ u = J2 

b k (p k satisfying b 2 k \\! 2 < +oo > , 

fc=l k=l ) 
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where d(x) = dist(x,dfl), and {Xki^k} is the Dirichlet spectral de- 
composition of —A in Q as above, with an orthonormal basis of 
L 2 (n). 

Furthermore, Vo(^) equipped with the norm 

IMk(n) = |||w|llri/a(n) + J ~Jj ( L9 ) 

is a Banach space. 

The fact that d~ l ' 2 u G L 2 (VL) if u is the trace of a function in 
Hq L (C) follows from a trace boundary Hardy inequality, originally due 
to Nekvinda [21] ; see Lemma 12.61 in next section for a simple proof. 
Thus, in next section we need to consider the operator A1/2 defined 
as in (II. 2p but now mapping A1/2 ■ Vo(fi) — > Vq(CI), where Vq(^) is 
the dual space of Vo(fi). For u = Y^k=ibk<Pk G Vo(f)), we will have 

= i ^fcV^fc) — X]fe=i ^fc^fe V^fc Moreover, there will be a unique 
harmonic extension v G Hq L (C) in C of w, and it is given by the ex- 
pression 

00 

v(x, y) = ^ h<Pk(x) exp(-\l /2 y) for all (x, y) G C. 
k=i 

Thus, the operator Ai/ 2 : Vo(O) — > Vq(^) is given by the Dirichlet- 
Neumann map 

o 00 

A I 7 xl/2 

Ai/ 2 « := t^j |nx{o}= 2^ bkA k Vk- 

k=l 

Note that A1/2 o A1/2 is equal to —A in Q with zero Dirichlet boundary 
value on dfl. More precisely, we will have that the inverse Bi/ 2 = A~, 2 
— which maps Vq(O) into itself, and also L 2 (Q) into itself — is the 
unique square root of the inverse Laplacian (— A) -1 in Q with zero 
Dirichlet boundary values on dfl; see next section for details. 

To establish the regularity of weak solutions to (II. ip obtained by the 
previous minimization technique, we establish the following results of 
Calderon-Zygmund and of Schauder type for the linear problem 

A 1/2 u = g(x) inn 

u = on oil; v ' 

see Proposition 13.11 for more details. 

Theorem 1.9. Let u G Vo(^) be a weak solution of (11.101) . where 
g G Vq{VL) and Q is a C 2 ' a bounded domain in W 1 , for some < a < 1. 
Ifge L 2 (Q), then u G H^(fl). 
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If ge HliVL), then u G # 2 (fi) D H%(Q). 
Ifge L°°(n), then u G C a (tt). 
Ifge C a (n^and g\ m = 0, then u G C^H). 
Ifge C^iU) and g\ dn = 0, then u G C 2 ' a (U). 

In this paper we will give full — and rather simple — proofs of these 
regularity results, specially since we could only find references for some 
of them and, besides, in close statements to ours but not precisely ours. 
Our proof of Theorem 11.91 uses the extension problem in Q x (0, oo) re- 
lated to fll.lOp . and transforms it to a problem with zero Dirichlet 
boundary in Q x {0} by using an auxiliary function introduced in [5]. 
Then, by making certain reflections and using classical interior regu- 
larity theory for the Laplacian, we prove Holder regularity for u and 
its derivatives. 

To apply the previous Holder regularity linear results to our nonlin- 
ear problem (11. ip . we first need to prove that g := f(u) is bounded, i.e., 
u is bounded. We will see that boundedness of weak solutions holds 
for subcritical and critical nonlinearities; we establish this result in sec- 
tion [5j We will follow the Brezis-Kato approach bootstrap method. In 
this way, we establish the following (see Theorem 15.21) . 

Theorem 1.10. Assume that go is a Caratheodory function in Q x K. 
satisfying 

\g (x, s)\ < C(l + \s\ p ) for all (x, s) G Q x R, 

for some constant C , 1 < p < ^j- if n > 2, or 1 < p < oo if n = 1, 
where Q is a smooth bounded domain in M. n . Let u G Vo(fi) be a weak 
solution of 

\A 1/2 u = g (x,u) inVt, 
1 u = on <9f2. 

Then, u G L°°(Q). 

The paper is organized as follows. In section [21 we study the ap- 
propriate function spaces Hi L {C) and Vo(^), and we give the proof 
of Proposition 11.81 and other related results. The regularity results of 
Theorem 11.91 can be founded in section [31 Maximum principles, Hopf's 
lemma, and the maximum principle in "small" domains of Proposition 
11.71 are proved in section [U The complete proof of Theorem II .11 is given 
in section [5] by studying the minimization problem and applying the 
previous results on regularity and maximum principles. We prove The- 
orem [TJTO] also in section [51 while Theorems 11.31 and II .51 are established 
in section [61 and Theorem 11.61 in section [71 
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2. Preliminaries: function spaces and the operator Ay% 

In this section we collect preliminary facts for future reference. First 
of all, let us set the standard notations to be used in the paper. We 
denote the upper half-space in IR n+1 by 

= {z = (x, y) = (an, ■ • • , x n , y) e | y > 0}. 

Denote by H S (U) = W S ' 2 (U) the Sobolev space in a domain U of M" 
or of M+ +1 . Letting U C E n and s > 0, iP(£/) is a Banach space with 
the norm 

f f f \u(x)- U (x)\ 2 J J f . . N|2j 

= ( J u J u ^—^^ dxdx + dx j ■ 

Let f2 be a bounded smooth domain in M n . Denote the half-cylinder 
with base Q by 

C = tt x (0,oo) 

and its lateral boundary by 

d L C = dVlx [0,oo). 

To treat the nonlocal problem (11. II) . we will study a corresponding 
extension problem in one more dimension, which allows us to investi- 
gate (11.11) by studying a local problem via classical nonlinear variational 
methods. We consider the Sobolev space of functions in H 1 ^) whose 
traces vanish on d^C: 

Hl L (C) = {v e H l {C) | v = a.e. on d L C }, (2.1) 

equipped with the norm 

\ 1/2 

\Vv\ 2 dxdy) . (2.2) 



We denote by tr^ the trace operator on fix {0} for functions in Hq L (C): 

tr n i> :=«(-,()), for veHl L (C). 

We have that ti n v G H 1 ' 2 ^), since it is well known that traces of H 
functions are H 1 ^ 2 functions on the boundary. 

Recall the well known spectral theory of the Laplacian —A in a 
smooth bounded domain Q with zero Dirichlet boundary values. We 
repeat each eigenvalue of —A in Q with zero Dirichlet boundary con- 
ditions according to its (finite) multiplicity: 

< Ai < A 2 < • • • < Afc < • • • — ► oo, as k — > oo, 
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and we denote by ipk G Hq(Q) an eigenfunction corresponding to A& 
for k — 1, 2, • • • . Namely, 

-A(p k = \ k (p k in tt, « 2 3 < 

</? fe = on f2. ^ ' ' 

We can take them to form an orthonormal basis {(pk} of L 2 (il), in 
particular, 

ip 2 k dx = 1, 



/ 

Jn 



in 

i2i 



and to belong to C (Q) by regularity theory. 

Now we can state the main results which we prove in this section. 

Proposition 2.1. Let Vo(Q) be the space of all traces on Q x {0} of 
functions in Hq L (C). Then, we have 

Vo(fi) := {u = tr n v\veHl L (C)} 

= \ueH 1 ' 2 {U)\ [ ^l d x <+oc\ 
L Jn d(x) J 

oo oo ^ 

u G L 2 {Q) I u = ^^bk^Pk satisfying ^^b^X^ 2 < +oo > , 

k=l k=l ) 

where d(x) = dist(x, dQ), and {Xk,<Pk} is the Dirichlet spectral de- 
composition of —A in Q as above, with {<fk} an orthonormal basis of 
L 2 (VL). 

Furthermore, Vq{VL) equipped with the norm 

,2 ^ 1/2 



\u 



u 



llvo(n) = { 11^11^/2(0) + / -J \ ( 2 - 4 ) 



is a Banach space. 



Proposition 2.2. If u G Vq{VL), then there exists a unique harmonic 
extension v in C of u such that v G Hq L (C). In particular, if the 
expansion of u is written by u = Y^iT=i bkfk G Vo(Q), then 

oo 

v ( x , V) = XI h ^k{x) exp(-\l /2 y) for all (x, y) G C, 

k=l 

where {\k, <fi>k} is the Dirichlet spectral decomposition of —A in Q as 
above, with {(pk} an orthonormal basis of L 2 (Q). 
The operator Ay 2 '■ Vo(fi) — > Vq(J1) is given by 

M/2U := ^~ 
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where Vg(fi) is the dual space ofVo(£l). We have that 

oo 

Ai/ 2 u = y^b k \]/ 2 Lp k , 



k=l 



and that A\/ 2 ° Ai/ 2 (when A 1 / 2 is acting, for instance, on smooth 
functions with compact support in Q) is equal to —A in Q with zero 
Dirichlet boundary values on dQ. More precisely, the inverse Bi/ 2 '■ = 
Ay 2 is the unique positive square root of the inverse Laplacian (— A) -1 
in Q with zero Dirichlet boundary values on dQ. 

The proofs of these two propositions need the development of sev- 
eral tools. First let us give some properties of the space Hq L (C). De- 
note by V 1 ' 2 (WY*~ ) the closure of the set of smooth functions com- 



pactly supported in IR™ +1 with respect to the norm of ||HI:di,2( R «+1) = 
( /r™+ 1 I Vw| 2 dxdy) . We recall the well known Sobolev trace in- 



equality that for w G V h2/mn+1 



+ J' 

(n-l)/2n ,t , Nl/2 



\ in— L) An /I \ 

|«>(ar,0)| an/(n - 1 W) <C[ / \Ww(x,y)\ 2 dxdy) 



where C depends only on n. 
Denote for n > 2, 



(2.5) 



2« = -*L and 2« -1 = ^ + 1. 
n — 1 n — 1 

We say that p is subcritical if 1 < p < 2* — 1 = for n > 2, and 
1 < p < 00 for n = 1. We also say that p is critical if p = 2^ — 1 = ^j- 
for n > 2, and that p is supercritical if p > 2$ — 1 = ^j- for n > 2. 
Lions 1 19 1 showed that 



f n +i \Vw(x,y)\ 2 dxdy 

is achieved. Escobar [12] prove that the extremal functions have all the 
form 

£ (n-l)/2 

C4(x, w) = 77 71 — 7, (2.7) 

where x$ G M n and e > are arbitrary. In addition, the best constant 
is 

(n — DaW 
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where a n denotes the volume of n-dimensional sphere S™ C M n+1 . 

The Sobolev trace inequality leads directly to the next three lemmas. 
For v G Hq L (C), its extension by zero in \C can be approximated 

by functions compactly supported in IR™ +1 . Thus the Sobolev trace 
inequality (I2.5P leads to: 

Lemma 2.3. Let n > 2 and 2" = -^r. Then there exists a constant C , 

— n— 1 ' 

defending only on n, such that, for all v G Hq L (C), 



By Holder's inequality, since Q is bounded, the above lemma leads 
to: 

Lemma 2.4. Let 1 < q < 2^ for n > 2. Then, we have that for all 



where C depends only on n,q, and the measure ofQ. Moreover, ( 12.91) 
also holds for l<q<ooifn=l. 

This lemma states that tr n (iJ 1 L (C)) C L«(0), where 1 < q < 2 fl 
for n > 2 and 1 < q < oo for n = 1, (see the proof of Lemma 12.51 
for the case n = 1). In addition, we also have the following compact 
embedding. 

Lemma 2.5. Let 1 < q < 2 tt = ^ for n > 2 and 1 < g < oo for 
n = 1. T/ien tr^^o l(^)) ^ s compactly embedded in L q (Q). 

Proof. It is well known that tr n ( J ff 1 L (C)) C H 1 ' 2 ^) and that H 1 / 2 ^) 
CC L 9 (fi) when 1 < q < 2 s = ^ for n > 2 and 1 < g < oo for n = 1. 
Here CC denotes the compact embedding. This completes the proof 
of the lemma. However, if one wants to avoid the use of the fractional 
Sobolev space if 1/,2 (f2), the following is an alternative simple proof. 

Considering the restriction of functions in C to Q x (0, 1), it suffices 
to show that the embedding is compact with C replaced by Q x (0, 1). 
To prove this, let v m G H^ L (tt x (0,1)) := {v G H\Q x (0,1)) | 
v = a.e. on <9f2 x (0, 1)} such that v m — weakly in L (f2 x (0, 1)), 
as m —>■ oo. We may assume by the classical Rellich's theorem in 
Q x (0, 1) that v m — > strongly in L 2 (Q x (0, 1)), as m — > oo. We 
introduce the function w m = (1 — y)v m . It is clear that 




(2.8) 





(2.9) 



w, 



'm|f2x{0} — V, 



rn ■ 



W 



'm\Qx{l} 



0. 
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By direct computations we have 

\v m (x,0)\ 2 dx = / \w m (x,0)\ 2 dx = - / d y (w 2 n (x, y)) dxdy 
n Jn Jo Jn 

- 2 {f o J w m( x >y) dxd y) {J Q J \^ w ™( x iv)\ 2 dxd v) ■ 

Therefore, since w m = (1 — y)v m is bounded in H l (Q x (0,1)) and 
w m — > strongly in L 2 (Q x (0, 1)), we find that, as m — > oo, 

v m (x, 0) — > strongly in L 2 (Q) and hence also in L 1 (f2). 

On the other hand, since q is subcritical, the following interpolation 
inequality, 

IM-,o)|U <(n) < \\v m {-M\W)\MM\ l J m 

for some < 6 < 1 completes the proof since we already know that v m 
converges strongly to zero in L 1 (f2). □ 

We also need to establish a trace boundary Hardy inequality, which 
already appeared in a work of Nekvinda [21] . 



Lemma 2.6. We have that 

tm(Hl L (C)) c H l l 2 (Sl) 
is a continuous injection. In addition, for every v G Hq L (C), 

dx<C I \Vv(x,y)\ 2 dxdy, (2.10) 
d[x) Jc 

where d(x) = dist(x, <9f2) and the constant C depends only on Q. 

Proof. The first statement is clear since the traces of H l (C) functions 
belong to H l l 2 (dC). Regarding the second statement, we prove it in 
two steps. 

Step 1. Assume first that n — 1 and Q = (0, 1). For < xq < 1/2, 
consider the segment from (0, xq) to (xq, 0) in C = (0, 1) x (0, oo). We 
have 

v(x Q ,0) = v(t,x - t) |^ = / (d x v - d y v)(t,x -t)dt. 

Jo 

Then 

\v(x ,0)\ 2 <x 2\Vv(t,x -t)\ 2 dt. 
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Dividing this inequality by xo and integrating in xq over (0, 1/2), and 
making the change of variables x = t, y = xq — t, we deduce 

1/2 I ( n)\2 rl/2 rl/2 , 

1 1 °' n dx <2 dx dy\Vv\ 2 <2 / \Vv\ 2 dxdy. 



x o Jo Jo 

Doing the same on (1/2,1), this establishes inequality f)2.10p of the 
lemma. 

Step 2. In the general case, after straightening a piece of the boundary 
dfl and rescaling the new variables, we can consider the inequality in 
a domain D = {x = (x',x n ) \ \x'\ < 1,0 < x n < 1/2} and assume 
that v = on {x n — 0, \x'\ < 1} x (0, oo), since the flatting procedure 
possesses equivalent norms. By the argument in Step 1 above, we have 



x n Jo JO 



dx n < C I / |Vf | dx n dy, 



for all x' with < 1. From this, integrating in x' we have 



I 



D x n J D Jo X n 

<C / \Vv\ 2 dxdy. 

JDx(0,od) 

Since after flattening of dfl, x n is comparable to d(x) = dist(x, dQ), 
this is the desired inequality (I2.10p . □ 

Recall that the fractional Sobolev space H l ^{ft) is a Banach space 
with the norm 

\\u\\ 2 m /2 l ^ = [ [ ^\ X ^ zt^tt dxdx+ [ \u(x)\ 2 dx. (2.11) 
(> JnJn \x-x\ n+1 J n 

1 /2 

Note that the closure H (fl) of smooth functions with compact sup- 
port, C c °°(fi), in H l / 2 {Vt) is all the space H l / 2 {Vt), (see Theorem 11.1 in 
[20]). That is, C£°(f2) is dense in f/" 1//2 (f2). However, in contrast with 
this, the trace in Q of functions in Hq L (C) "vanish" on dQ in the sense 
given by f[2TT0]) . 

Recall that we have denoted by Vq(Q) the space of traces on Q x {0} 
of functions in Hq L (C): 

Vo(fi) := {u = ti n v | v E Hl L {C)} C F 1/2 (fi), (2.12) 

endowed with the norm (12 .4p in Proposition 12.11 The dual space of 
Vo(fi) is denoted by Vg(f2), equipped with the norm 

IMIv*(n) = sup{(u,#) I u E Vo(fi), ||ti||vb(n) < !}• 
Next we give the first characterization of the space Vq(Q): 
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Lemma 2.7. Let Vq(£1) be the space of traces on Q x {0} of functions 
in Hq L (C), as in ( I2.12p . Then, we have 

Vo(fi) = {ue tf 1/2 (fi) I f ^-dx < +oo) , 
where d(x) = dist(x, d£i). 

Proof. The inclusion C follows Lemma [2.61 Next we show the other 
inclusion. Let u G if 1/,2 (f2) satisfy J Q u 2 /d < oo. Let u be the extension 
of u in all of W 1 assigning u = in R n \ Q. The quantity 

\\u\\ 2 Hl/ 2 m = I f ^[^ dxdx+ I \u{x)\ 2 dx 

can be bounded — using u = in IR n \ Q — by a constant times 

\ u \\ 2 m/ 2 (Q,) + I" ~ ~ dx 



d(x) 



that we assume to be finite. Hence, u £ if 1/,2 (IR n ) and thus u is the 
trace in R n = dRl +1 of a function v G H 1 (Rl +1 ). 

Next, we use a partition of the unity, and local bi-Lipschitz maps 
(defined below) sending M" +1 into f2 x [0, oo) = C being the identity 
on Q x {0} and mapping R n \ tt = (9M" +1 ) \ tt into d£l x [0, oo). By 
composing these maps with the function (cutted off with the partition 
of unity) v, we obtain an Hq L (C) function with u as trace on Q x {0}, 
as desired. 

Finally we give a concrete expression for one such bi-Lipschitz maps. 
First, consider the one dimensional case Q = (0, oo). Then simply take 
the bi-Lipschitz map 

(x, y) G (0, oo) x (0, oo) = Q x (0, oo) 

. x 2 — y 2 2xy . 



a/x 2 + y 2 \/ x 2 + y 2 



x (0, oo) 



whose Jacobian can be checked to be identically 2. In the general case, 
we can flatten the boundary dfl and use locally the previous map. □ 

Next we consider, for a given function u G Vq(Q), the minimizing 
problem: 

inf | J | Vtf dxdy \ v G ifJ L (C), 0) = u in ft j . (2.13) 

By the definition of Vq(Q), the set of functions v where we minimize is 
non empty. By lower weak semi-continuity and by Lemma 12.51 we see 
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that there exists a minimizer v. We will prove next that this minimizer 
v is unique. We call v a weak solution of the problem 

Av = in C, 

v = ond L C, (2.14) 
v = u on Vt x {0}. 

That is, we have 

Lemma 2.8. For u G Vo(0) 7 there exists a unique minimizer v of 
(I2.13p . The function v G Hq L (C) is the harmonic extension of u (in 
the weak sense) to C vanishing on dijZ. 

Proof. By the definition of Vq(Q), we have that, for every u G Vo(Q), 
there exists at least one w G Hi L {C) such that tr^(w) = u. Then the 
standard minimization argument gives (using lower semi-continuity and 
Lemma |2~5]) the existence of a minimizer. The uniqueness of minimizer 
follows automatically from the identity of the parallelogram used for 
two possible minimizers v\ and t>2, 

o < J(^) = \JM + \jm - J(^) < o, 

where J(v) = f c \ Vf | 2 dxdy, which leads to v± = v 2 - □ 

By Lemma [2.81 there exists a unique function v G Hq L (C) which is 
the harmonic extension of u in C vanishing on 8lC, and that we denote 
by 

v := h-ext(u). 

It is easy to see that for every r\ G C°°(C) fl H 1 ^) and 77 = on d^C, 

AV vV r] dxdy = [ ^-r/dx. (2.15) 

By Lemma [2761 there exists a constant C such that for every u G Vo(fi), 

||w||v (n) < C||-h-ext(?x)|| H i i{c) . (2.16) 

Next, note that the h-ext operator is bijective from Vq(O) to the sub- 
space Tt of Hq L {C) formed by all harmonic functions in Hq L (C). Since 
both Vo(^) and 7i are Banach spaces, the open mapping theorem gives 
that we also have the reverse inequality to (12.161) . i.e., there exists a 
constant C such that 

\\h-ext(u)\\ HlL{c) <C\\u\\ Vo(nh (2.17) 

for all u G Vo(O). From this we deduce the following. Given a smooth 
£ G Vo(O), consider the h-ext (£) and call it 77. Now, we use (12 . 15[) and 
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f l2~T7j) (for u and £) to obtain / n dx 
is, |nG Vq(H) and there is the bound: 



< C\\u\ 



Vo(Q) 



That 



— — h-ext(-u) 

<7Z/ 



< C||u| 



v *(n) 



v (n)- 



Hence we have 

Lemma 2.9. JTie operator Ai/ 2 : Vq(^) — >■ Vg(fi) defined by 



dv 



(2.18) 



nx{o} 



where v = h-ext(u) G Hq L (C) is the harmonic extension of u in C 
vanishing on 8lC, is linear and bounded from Vo(fi) to Vq(Q). 

We now give the spectral representation of Ai/ 2 and the correspond- 
ing structure of the space Vq{VL). 

Lemma 2.10. (i) Let {<p k } be an orthonormal basis of L 2 (Q) forming 
a spectral decomposition of —A in Q with Dirichlet boundary conditions 
as in (12. 3p . with {X k } the corresponding Dirichlet eigenvalues of —A 
in Q. Then, we have 



I k=l k=l J 

(ii) Let u G Vq(Q). Then we have, if u = J2T=i bktyk, 

oo 

A 1/2 u = $3& fc Ai /2 V fc eV *(fi). 

k=l 

Proof. Let u G Vo(fi), which is contained in L 2 (Q). Let its expansion 
be written by u(x) = YlT=i bk^kix)- Consider the function 



v(x,y) = ^2hip k (x)exp(-\l /2 y), 



(2.19) 



k=l 



which is clearly smooth for y > 0. Observe that v(x, 0) = u(x) in Q 
and, for y > 0, 



El /2 1/2 
b k {-X k ip k (x) exp(-A fc 7 y) + X k <f k (x)exp(-X k ' y)} = 



k=l 



Thus, v is a harmonic extension of w. We will have that v = h-ext(w), 
by uniqueness, once we find the condition on {b k } for v to belong to 
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Hq L {C). But such condition is simple. Using ( 12. 19ft and that {<fk} are 
eigenfunctions of —A and orthonormal in L 2 (fl), we have 

/*oo i* 

\Vv\ 2 dxdy = I {\V x v\ 2 + \d y v\ 2 } dxdy 
Jo Jn 

OO nQQ 

= 2j2b 2 k Xk / exp(-2\l /2 y)dy 

JO 



fc=l 
OO 



k=l k=l 

This means that v G Hq L (C) if and only if Ylk^i^k^k 2 < 00 ■ There- 
fore, this condition on is equivalent tone Vo(Q). 

Assertion (ii) follows from the direct computation of — ^ \ y= o using 

In functional analysis, the classical spectral decomposition holds for 
self-adjoint compact operators, such as the Dirichlet inverse Laplacian 
(—A) -1 : L 2 (Q) — > L 2 (Q). This is the reason why we now define, with 
the aid of the Lax-Milgram theorem, a compact operator B\i 2 which 
will be the inverse of Ai/ 2 . 

Definition 2.11. Define the operator Bi/ 2 : Vq(O) — > Vo(fi), by g t— > 
tr^f , where v is found by solving the problem: 

Av = in C, 

u = on d L C, (2.20) 
g = </(x) onfix{0}, 

as we indicate next. 

We say that v is a weak solution of (12.201) whenever v G Hq L (C) and 

VvV£dxdy = (g,£(;0)) (2.21) 

for all £ G Hq L (C). We see that there exists a unique weak solution 
of (12.201) by the Lax-Milgram theorem, via studying the corresponding 
functional in Hq l (C): 

I{v) = ^JjVv\ 2 dxdy- (g,v{-,0)), 

where g G Vg(f2) is given. Observe that the operator B\/ 2 is clearly the 
inverse of the operator Ai/ 2 . 

On the other hand, let us compute By 2 o By^^^y Here note that 
since V (fi) C L 2 (tt), we have L 2 (fi) C V£(fi). For a given # G L 2 (fi), 
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let ip G H^ip) n H 2 (fl) be the solution of Poisson's problem for the 
Laplacian 

—A<p = g in Q, 
ip = on dQ. 

Since Hq(Q) C Vo(fi) (for instance, by Lemma 12.101) . there is a unique 
harmonic extension ip G Hq L (C) of ip in C such that 

= in C, 
•0 = on diC, 
ip = (p on f2 x {0}. 

Moreover, ip(x,y) := ip(x,y) — <p(x) solves 

Atp = —g(x) in C, 

on d^C, 
on x {0}. 

Considering the odd reflection ip od of ^ across f2 x {0}, and the function 

9od(x,y) = 





g{x), y>0, 

-g(x), y < o, 

we have 

-g a d in fi x I, 
on dVl x R. 

Therefore, since g Q d G L 2 (Qx (— 2, 2)), we deduce G H 2 (Qx (— 1, 1)) 
and hence ^ G /f 2 (fi x (0, 1)). We deduce, by the smoothness of the 
harmonic function ip for y > and by its exponential decay in y — see 

( ESSD — , that v e ^o,l( c ) n # 2 (C). 

It follows that —dyip G Hq L (C) solves 

A(-d y ip) = in C, 
— = on d L C, 



-dyip) = dy y ip = - A x ip = - A(p = g on Q x {0}. 



and 

Since V Q (fi) C L 2 (fi), we have that # G L 2 (fi) = L 2 (ft)* C V *(fi), 
and we deduce that the solution v G Vo(fi) of (I2.20p is v = —d y ip, 
because of the uniqueness of Hq l {C) solution of f)2.20p . In particular, 
B\/29 = ^(-,0) = — d y ip(-, 0). On the other hand, since ip G Hq L (C) 
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solves 

' Atp = in C, 

ip = on 8lC, 

§£ = -d y il;(; 0) = v(; 0) = S 1/2 p on Q X {0}, 

we conclude that 

(Si /2 o S 1/2 )(/ = £? 1/2 <, 0) = 0) = y> = (-A)"V 

Summarizing the above argument, we have: 

Proposition 2.12. B 1/2 o B 1/2 \ L2(n) = (-A)" 1 : L 2 (Vt) -> L 2 (fi) ; w/jere 
(—A) -1 zs i/ie inverse Laplacian in Q with zero Dirichlet boundary 
conditions. 

Note that B\i 2 '■ L 2 (Q) — > L 2 (f2) is a self-adjoint operator. In fact, 
since for t>i,t> 2 G Hq L (C), 

(v 2 Avi — ViAv 2 ) dxdy = [ (v 2 ^^- — Vi^-^-) dx, 



dv dv 



we sec 



and 



/ Bi /2 g 2 - g x dx = I B 1/2 g 1 -g 2 dx 
Jn Jn 



v 2 (x, 0)A l/2 v 1 (x, 0) dx — / 1)1(0;, 0)A 1/2 v 2 (x, 0) Gfo. 

On the other hand, by using (12.211) with £ = v and Lemma I2.5[ we ob- 
tain that Bu2 is a positive compact operator in L 2 (Q). Hence by the 
spectral theory for self-adjoint compact operators, we have that all the 
eigenvalues of B\j 2 are real, positive, and that there are corresponding 
eigenfunctions which make up an orthonormal basis of L 2 (Q). Further- 
more, such basis and eigenvalues are explicit in terms of those of the 
Laplacian with Dirichlet boundary conditions, since (— A) -1 has Bu 2 
as unique, positive and self-adjoint square root, by Proposition 12.121 
Summarizing: 

Proposition 2.13. Let be an orthonormal basis of L 2 (Q) form- 
ing a spectral decomposition of —A in Q with Dirichlet boundary con- 
ditions, as in (12.31) . with {\k} the corresponding Dirichlet eigenvalues 
of —A in Q. Then, for all k > 1, 

Ai/ 2 <fik = K /2( Pk m Vt, / 2 
ip k = on dQ. 

In particular, {(pk} is also a basis formed by the eigenfunctions of A 1 / 2 , 
with eigenvalues {A^ 2 }. 
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Proof of Proposition |2. li lt follows from Lemma and Lemma ^.lOl 

□ 

Proof of Proposition 12.21 It follows from Lemma 12.81 Lemma 12.101 
and its proof, and Propositions 12.121 and 12.131 □ 



3. Regularity of solutions 

In this section we study the regularity of weak solutions for linear 
and nonlinear problems involving Ai/ 2 - First we consider the linear 
problem 

f A l/2 u = g(x) in SI, . g ^ 

I ii = fl on dfl, ^ ' ' 

where g G Vq(Q) and Q is a smooth bounded domain in M. n . By the 
construction of the previous section, the precise meaning of (13.11) is 
that u = tiQV, where the function v G Hq L (C) with v(-, 0) = iiG Vo(^) 
satisfies 

Av = in C, 

v = on d L C, (3.2) 
% = g{x) onfix{0}. 

We will say then that v is a weak solution of (13. 2ft and that u is a weak 
solution of (13.11) . 

Most of this section contains the proof of the following analogues 
of the iy 2 ' p -estimates of Calderon-Zygmund and of the Schauder esti- 
mates. 

Proposition 3.1. Let a G (0, 1), Vt be a C 2)OL bounded domain ofW 1 , 
g G Vq(Q), v G Hq L (C) be the weak solution of (13.21) . and u = tr^v be 
the weak solution of (13. ip . Then, 

(i) Ifge L 2 (Q), then u G H^(D). 

(ii) Ifge Hl(tt), then u G H 2 (tt) n H%(Q). 

(iii) If g G L°°(Q), then v G W^> q (tt x (0,R))Jor all R > and 
1 < q < oo. In particular, v G C a (C) and u G C a (f2). 

(iv) Ifge C a (jT) and g\ 9n = 0, then v G C x ' a (C) and u G C^H). 

(v) Ifge C l ' a (yi) and g\ dn = 0, then v G C 2 > a (C) and u G C 2 ' a (U). 

As a consequence, we deduce the regularity of bounded weak solu- 
tions to the nonlinear problem 



Ai/ 2 u = f(u) in ft, 
u = on dfl. 



(3.3) 
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As before, the precise meaning for ( 13. 3ft is that v G Hq L (C), v(-, 0) = u, 
and v is a weak solution of 

Av = in C, 

u = on d L C, (3.4) 

£ = /(«(-, 0)) onfix{0}. 

Here the weak solution u is assumed to be bounded. Regularity results 
for weak solutions not assumed a priori to be bounded, of subcritical 
and critical problems will be proved in section 

By Cq(Q) we denote the space of continuous functions in Q vanishing 
on the boundary dQ. In the following result note that /(0) = is 
required to have C 1 (fi) regularity of solutions of (13.31) . 

Proposition 3.2. Let a G (0, 1), Q be a C 2 ' a bounded domain ofW 11 , 
and f be a C 1,a function such that /(0) =0. // u G L°°(Q) is a weak 
solution of ([E3D, and thus v G H% L (C) H L°°(C) zs a u>ea£; solution of 
(I3T4D . toen u G C 2 - Q (n) n C (n). /n addition, v G C 2 ' a (C). 

Proof. By (hi) of Proposition 13.11 we have that u G C a (Q). Next, by 
(iv) of Proposition l3.1l and since on <9f2x {0}, g := /(f (•, 0)) = /(0) = 0, 
we have u G C 1 '"^). Finally, v G C 2 ' a (C) and u G C 2 ' a (H) from (v) of 
Proposition 13.11 since g = f(u) vanishes on dfl and it is of class C 1,a , 
since both / and u are C 1,a . 

□ 

Proof of Proposition [3711 (i) and (ii). Both statements follow im- 
mediately from Propositions 12.11 and 12.21 Simply use that {(fk} is an 
orthonormal basis of L 2 (f2) and that {(fk/ ^]/ 2 } is an orthonormal basis 
of Hq(Q). For part (ii), note that if A\jiu — g G Hq(Q), then we have 

AueL 2 (n). 

(iii) Let v be a weak solution of (13. 2p . We proceed with a useful 
method, introduced by Cabre and Sola Morales in [5], which consists 
of using the auxiliary function 

ry 

w(x,y) — / v(x,t)dt for (x, y)EC. (3.5) 
Jo 

Since (Aw) y = in C, we have that Aw is independent of y. Hence 
we can compute it on {y = 0}. On {y = 0}, since w = 0, we have 
Aw = w yy = Vy. Thus w is a solution of the Dirichlet problem 



f -Aw(x,y) = g(x) in C, 
\ w = on dC. 



(3.6) 
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We extend w to the whole cylinder H x R by odd reflection: 

w(x, y) for y > 0, 
— w(x, —y) for y < 0. 



y) = 

Moreover, we put 

9od(x, y) 

Then we obtain 



g(x) for y > 0, 
—g(x) for 1/ < 0. 



-Au; od = g od in x R 
Wod = on <9f2 x 



(3.7) 



9h 




Since g 0(i G L q (Q x (— 2f?, 2f?)) for all R > and 1 < g < oo, regularity 
for the Dirichlet problem (13.71) gives w od G W /2 ' 9 (fi x (-R,R)) for all 
f? > and 1 < g < oo. In particular, w G C 1,a (C). Therefore, 
v = w y G C a (C) and u G C a (H). 

(iv) Choose a smooth domain if such that Q C H, and let 

g in f2, 
in if \ a 

We have that gn G C a (ff), since g|dn= 0, by assumption. Consider 
the weak solution of 

At># = in if x (0, oo), 
on dH x [0,oo), 

= ^(ac) on if x {0}. 

Consider also the auxiliary function 

rv 

wh{x,v) = / vn(x,t)dt in if x [0, oo), 
Jo 

which solves problem (13.61) with f2 and g replaced by if and 

Using boundary regularity theory (but away from the corners of if x 
[0, oo)) for this Dirichlet problem, we see that % is C 2 ' a (H x (0, oo)) 
(again, here we do not claim regularity at the corners dH x {0}). Thus, 
wh G C 2 ' a {C) (here instead we include the corners dVl x {0} of C). 

Consider the difference ip = wh — w in C, where w is defined by (13. 5p . 
It is clear that 

' Ay? = in C, 
ip = w H on d L C, 
<p = on Q x {0}. 

We extend ip to the whole cylinder O x 1 by odd reflection: 
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Moreover, we put 



Then we have 




A(p Q d = in fixi, 
Pod = w H ,od on dQ x R. 



w H (x,y) fory>0, 
—wh(x, —y) for y < 0. 



(3.8) 



Since G C 2 ' a (C), iujj = on dfl x {0}, and d yy WH = d y vn = 
—g H = — g = on dQ x {0}, we deduce that w H ,od G C 2,a (dtt x R). It 
follows from elliptic regularity for (13.81) that ip od G C 2 ' a (fi x R). Thus, 
ip G C 2 ' Q (C), u; G C 2 ' Q (C) and t> = d y w G C^C). _ 

(v) Choose a smooth bounded domain B such that C B. B 
could be the same as H in (ii), for instance a ball, but we change its 
name for notation clarity. Since g G C 1,a (n), there exists an extension 
Qb £ C 1,Q! (£>); see [15]. Consider the solution vb of 



As before, from interior boundary regularity for the Dirichlet prob- 
lem of the type (13.61) satisfied by wb, we obtain that wb G C 3,a (B x 
[0, oo)) since g B G C 1,a (B) (away from the corners dB x {0}). Thus, 
v B G C 2 ' a (B x [0,oo)). Thus, v B G C 2 ' a (C). Consider the difference 
ip = vb — v in C, where v is a weak solution of (13. 2p . We have that 
ip = vb — v satisfies 



We extend ip to the whole cylinder f2 x R now by even reflection: 




Consider the auxiliary function 





Atp = in C, 
ip = v B on d L C, 
f£ = on J] x {0}. 



ipev(x,y) = 



ip(x, y) for y > 0, 
?/>(x, — y) for y < 0. 



Moreover, we put 



v B ,ev(x,y) = 



v B (x,y) fory>0, 
v B (x,-y) fory<0. 
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Then, since ^ = on x {0}, we have 

f AV'et, = in fixl, 

\ i>ev = V B ,ev Oil 8Q X R. 



Since vb E C 2,a (C), —d y VB = 9b = 9 = on dQ x {0}, we deduce that 



In this section we establish several maximum principles for Ay 2 . We 
denote by Cq{VL) the space of continuous functions in Q vanishing on 
the boundary dQ. For convenience, we state the results for functions 
in Cq(Q) fl C 2 {VL) (a space contained in Hq(Q) C Vo(Q)), but this can 
be weakened. 

The first statement is the weak maximum principle. 

Lemma 4.1. Assume that u E C 2 (Q) satisfies 



where Q is a smooth bounded domain in R n and c > in fi. Then, 
u > in Q. 

Proof. Consider the extension v = h-ext(w). If we prove that v > 
in C, then u > in Q. Suppose by contradiction that v is negative 
somewhere in C. Then, since Av = in C and v = on d L C, we deduce 
that v is negative somewhere in Q x {0} and that infe v < is achieved 
at some point (x ,0) E Q x {0}. Thus, we have 



V B ,ev £ C 2 '°(C 

that ^ e C* 2 '' 




4. Maximum principles 




inf v = v(xq, 0) < 0. 



By Hopf 's lemma, 



v y (x o ,0) > 0. 



It follows 

dv 
dv 

Therefore, since c > 



v y (x , 0) = A 1/2 v(x , 0) < 0. 



^1/2^(^0, 0) + c(x )v(x , 0) < 0. 
This is a contradiction with the hypothesis Ai/ 2 u + c(x)u > 0. 
The next statement is the strong maximum principle for Ay 2 - 



□ 
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Lemma 4.2. Assume that u G C 2 (Q) satisfies 

A1/2U + c(x)u > in Q, 
u > in Q, 

u = on dQ, 

where Q is a smooth bounded domain in MJ 1 and c G L°°(Q). Then, 
either u > in Q, or u = in O. 

Proof. The proof is similar to that of Lemma 14.11 Consider t> = 
h-ext(n). We observe that v > in C. Suppose that v ^ but it = 
somewhere in Q. Then there exists a minimum point (xq, 0) G f2 x {0} 
of v where v(xq, 0) = 0. Then by Hopf's lemma we see that A 1 / 2 u(xo) = 
—v y (x ,0) < 0. This implies that Ai/ 2 u(x ) + c(x )u(x ) < 0, because 
of v(xo, 0) = u(xo) =0. □ 

Next we establish a Hopf lemma for Ax/2, following a proof from [9]. 

Lemma 4.3. Let f2 6e a bounded domain in W n and c G L°°(Q). 

(i) Assume that Q is smooth and that O^iiG C 2 (fi) satisfies 

A\nu + c(a;)u > in f2, 
it > in Q, 

u = on <9fi. 

Then, < on <9£7 ; where Uq is the unit outer normal to dil. 

(ii) Assume that P G <9fi and taat <9f2 is smooth in a neighborhood 
of P. LetO^ve C 2 (C) n L°°(C), waere C = SI x (0, 00), sate/?/ 

At> = in C, 

v > on <9^C, 

g + c(x)t;>0 onfix{0}. 

Ifv(P, 0) = ; £/&en ^j^ -* < 0, where u is the unit outer normal in M. n 

to on. 

Proof. We follow the proof given in [9]. Note that statement (i) is a 
particular case of (ii). Thus, we only need to prove (ii). 

Step 1. We shall first prove the lemma in the case c = 0. Without 
loss of generality we may assume that (P, 0) = Pi = (61, 0, • • • , 0) G 
dfl x {0}, 61 > and vq = (1, 0, • • • , 0). Hence we need to prove 

dv(Pi) <Q 

Since is smooth in a neighborhood of P, there is a half-ball in 
included in the domain C, such that Pi is the only point in the closed 
half-ball belonging also to d L C. Let P 2 G fi x {0} and r > be the 
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center and radius of such ball. Then we have P 2 = (b 2 , 0, ■ • • , 0) G 
Qx {0}. Denote 

B+{P 2 ) := {z = (x, y)\\z-P 2 \< \Pt -P 2 \=:r,y>0}cC, 
B+ /2 (P 2 ) := {z = (x, y)\\z-P 2 \< \P X -P 2 \/2,y> 0}, 

A = B+(P 2 )\B+ 2 (P 2 ). 

Recall that P x G dB+{P 2 ) n (dQ x {0}). 
Consider the function on A: 

<p(z) = exp(-A|z - P 2 \ 2 ) - exp(-A|Pi - P 2 | 2 ), 

with A > to be determined later. Note that 

Aip = exp(-A|z - P 2 \ 2 ) {A\ 2 \z - P 2 \ 2 - 2(n + 1)A} . 

We can choose A > large enough such that Aip > in A. 

On the other hand, by Lemma 14.21 we see that v > in A \ {Pi}. 
Hence, since if = on dB+(P 2 ) n {y > 0}, we can take e > such that 

v - ELp > on dA n {y > 0}. 

Since — A(v — E(p) > in A, and 

dv 

-d y (v-£ip) = — >0 ondA(l{y = 0}, 

(recall that c = 0) by the maximum principle as in Lemma 14.11 we 
obtain 

v — etp > in A. 

Thus, from v—eip = at P x we see that d Xl {v — e<p)(Pi) < 0. Therefore, 
d xl v(Pi) < ed Xl ip{P{) = -2A(6i -6 2 )e- A l Pl - P2 l 2 < 0. Thus we have the 
desired result. 

Step 2. In the case c ^ 0, we define the function w = v exp(— (3y) for 
some (3 > to be determined. From a direct calculation, we see that 

-Aw - 2f3d y w = (3 2 w > in C 

and, choosing f3 > \\c\\L°°(n), 

-d y w >[f3- c(x)]w > on Q x {0}. 

Now we can apply to w the same approach as in Step 1 , with A replaced 
by A + 2j3d y , and obtain the assertion. □ 

Finally, we establish a maximum principle for A x / 2 in domains of 
small measure. Note that in part (ii) of its statement, the hypothesis 
on small measure is made only on the base of Q of the cylinder C. 
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Proposition 4.4. (i) Assume that u G C 2 (Q) satisfies 

Ai/ 2 u + c(x)u > inQ, 
u = on dfl, 

where Q is a smooth bounded domain in IR n and c G L°°(Q). Then, 
there exists 5 > depending only on n and \\c~\\loo^), such that if 
|fi n {u < 0}| < 5, then w > in Q. 

(ii) Assume that Q is a bounded (not necessary smooth) domain of 
W 1 and c G L°°(Q). Let u G C 2 (C) n L°°(C) ; w/iere C = ft x (0, oo) ; 
satisfy 

Av = in C, 

f > on 
jg + c(:r)t;>0 onfix{0}. 

T/ien, t/iere exzsfo 5 > depending only on n and \\c~\\l°°(p), such that 
if\nn{v(-,0)<0}\<8thenv>0inC. 

Proof. For part (i) of the theorem, consider v = h-ext(w). We see 
that v satisfies the assumptions on part (ii) of the theorem. Hence, it 
is enough to prove part (ii). For this, let v~ = max{0, — v} > 0. Since 
v ~ = on dil x [0, oo), we see 

= / v~ Av dxdy = f v~*^-dx + f |Vt> ~\ 2 dxdy. 
Jc Jnx{o} uv J c 



Then, 

?! — 1^ rl/rrlii = — / u~ 

dv 



.... . , f _dv , 

vv I dxdy = — v — dx 



S]x{0} 



/ v cvdx= —c(v ) 2 dx 
JQx{0} Jq 

I c-{v-(-,0)) 2 dx 
Jnr\fv-<-,o)>o\ 



< 

'nn{i)-(-,o)>o} 

<|finK(.,o)>o}| 1 /™|| c -||^ (n) || w -(.,o)||i 2 „ /( „_ 1)(n) . 

Thus, extending v~ by outside C we obtain an i? 1 (R" +1 ) function 
and thus we have 



/ R n+i|Vu \ 2 dxdy j c \Vv~\ 2 dxdy 

\\ V ~('l 0)||^2„/(n-l)( K n) \\ V ~('l 0)||^2n/(n-l)(n) 

<|finK(-,0)>0}| 1 /"||c-|| L o O(f , ) , 



where So is the best constant of the Sobolev trace inequality in 

If \Q fl {v ~(-, 0) > 0}| is small enough, we arrive at a contradiction. □ 
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5. SUBCRITICAL CASE AND L°° ESTIMATE OF BREZIS-KATO TYPE 

In this section, we study the nonlinear problem (11.11) with f(u) = u p 
in the subcritical and critical cases. In the subcritical case we look for 
a function v(x, y) satisfying for and y G R + , 

Av = in C = fix(0,oo), 
v — on d L C = dQ x [0, oo), , , 

% = v p on n x {0}, ^ ^ 

v > in C, 

where v is the unit outer normal to C at Q x {0} and 1 < p < 2" — 1 
if n > 2, orl<p<ooifn = l. If t> is a solution of (15. ip . then 
v(x, 0) = u(x) is a solution of (11.11) with the nonlinearity f(u) = u p . 

In order to find a solution of (15 .ip as stated in Theorem 11.11 we 
consider the following minimization problem: 

/ = inf|y \Vv(x,y)\ 2 dxdy\veHl L (C),J \v(x, 0)\ p+1 dx = 

We show that Iq is achieved. 

Proposition 5.1. Assume that 1 < p < 2^ — 1 if n > 2 or 1 < p < oo 

if n = 1. T/ien Jo achieved in Hq l (C) by a nonnegative function v. 

Proof. First, there is a function v G L (C) such that 



/ |Vf (x, y)\ 2 dxdy < oo and / \v (x, 0)\ p+1 dx 
Jc Jn 



1. 



In fact, it suffices to take any C°° function with compact support in 
Q x [0, oo) and not identically zero on Q x {0}, and multiply it by 
an appropriate constant. Next we complete the proof by weak lower 
semi-continuity of the Dirichlet integral and by the compact embedding 
property in Lemma 12.51 Finally, note that \v\ > is a nonnegative 
minimizer if v is a minimizer. □ 

To establish the regularity of the minimizer just obtained, we prove 
an L°°-estimate of Brezis-Kato type by the technique of bootstrap for 
subcritical or critical nonlinear problems. Let go be a Caratheodory 
function in Q x R satisfying the growth condition 

\g (x,s)\ < C(l+ \s\ p ) for all (x, s) G SI x R, (5.2) 

where fl is a smooth domain in IR n , 1 < p < if n > 2, or 1 < p < oo 
if n = 1. We consider the problem 

Av = in C = flx (0,oo), 

v = on d L C = dQ x [0, oo), (5.3) 

& = go(;v) on nx{0}. 
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Theorem 5.2. Letv G Hq L (C) be a weak solution of ( 15. 3ft and assume 
the growth condition (15. 2\i for g , with 1 < p < ^j- if n > 2, 
\<p<ooifn = l. Then, v(-,0) G 



or 



Proof. The proof follows the one of Brezis-Kato for the Laplacian. 
First of all, let us rewrite the condition on go as 

\g Q (x,v)\ < a(x)(l + \v(x,0)\) 

with a function 

. = \g (x,v(x, 0))| 
1 J ' 1 + |v(ar, 0)| 

which satisfies 

< a < C(l + |u(-,0)| p - 1 ) G L n (fi), 

since v G H^ L (C), v{-,0) G L^(fi) and p - 1 < 
Denote 

= {(x,v) I \(x,y)\ < r and y > 0}. 

For /3 > and T > 1, let (p^x = vVj? G Hq L (C) with t>y = min{|t>|,T}. 
Denote 

D T = {(a:,y)eC| \v(x, y)\ < T}. 
By direct computation, we see 

J |V(ff^)| 2 dxdy = J nf\ Vt> \ 2 dxdy 

+ [ (2/5 + (3 2 )\v\ w \Vv\ 2 dxdy. 
Jd t 

Multiplying (15.3j) by y^T and integrating by parts, we obtain 
v^\Vv\ 2 dxdy + 2(3 \v\ w \Vv\ 2 dxdy 

JC J D T 

= / VvV(vv%?) dxdy = / g (x,v)vVj? dx 
Jc ' Jnx{o} 

< I a(x)(l + \v\fvfdx. 

JQx{0} 

Combining these facts, we have 

/ \V{vv^)\ 2 dxdy <C{(3 + 1) I a(x)(l + \v\ 2 )vf dx, 
Jc ' Jnx{o} 
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where C denotes different constants independent of T and of (5. By 
Lemma [23, we deduce 

/ f \ 2/2" r 

( / \vv^f dx) <C(p + l) / a(x)(l + \v\ 2 )vfdx. (5.4) 

Assume that |w(-,0)|^ +1 G L 2 (Cl) for some /3 > 0. Then we obtain 
that J nx { | H 2 fy ^ x an d Jsixfo} are bounded uniformly in T. In 

what follows, let C denote constants independent of T — but that may 
depend on (3 and ||t>(-, 0)^ +1 1| z, 2 (fi) • Given M > 0, we have 

a\v\ 2 v?f dx < Mq I \v\ 2 v?f dx + / a\v\ 2 v^ dx 

Qx{0} ' Jnx{0} J{a>M } 

< CM +[ I a n dx) Vn ( [ \vv^\ 2t dx" 2 ' 1 



'oil I a a x 1 1 ; \ vv t\ 

'{a>AI } ' Wfix{0} 



Wf7x{0} 7 



where e(M ) = (j^ a>M ^a n dx) 1 /™ — > as M — >■ 00. Note that we 
can deal with / nx r i clv^ dx in the analogue procedure. Therefore, we 
deduce from the last inequalities and (15.4)) . taking M large enough so 
that C(J3 + l)e(Mo) = |, that 

/ r \ 2/2" 

( / \w4fdx) <C(1 + M ). (5.5) 

Wftx{0} ' 

Thus letting T — >■ 00, since C is independent of T, we obtain that 
|u(-,0)|^ +1 G L 2tt (fi). This conclusion followed simply from assuming 
\v{-,0)f +1 G L 2 (Cl). 

Hence, by iterating (3q — 0, $ + 1 = + l)rz? if i > 1 in 

(15.51) . we conclude that v(-,0) G L 9 (f2) for all q < 00. Finally, the 
proof of part (iii) in Proposition 13.11 — which only uses g G L q (Cl) for 
all q < 00 and not g G L°°(Cl) — applied with g(x) = go(x, v(x, 0)), 
which satisfies \g\ < C(l + |w(-,0)| p ) G L g (Cl) for all g < 00, leads to 
v(-,0) G C a (n) C L°°(Cl). □ 

Proof of Theorem 11.11 Proposition 15.11 gives the existence of a 
weak nonnegative solution v to (15. ip after multiplying the nonnegative 
minimizer of Iq by a constant to take care of the Lagrange multiplier. 
Then, Theorem 15.21 gives that f(-,0) G L°°(Cl). Next, Proposition 13.21 
gives that u G C 2,a (Cl), since f(s) = \s\ p is a function for some 
a G (0, 1). Finally, the strong maximum principle, Lemma |4.2[ leads 
to u > in Cl. □ 
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6. A PRIORI ESTIMATES FOR POSITIVE SOLUTIONS 

In this section we prove Theorem 11.31 Namely, we establish a priori 
estimates of Gidas-Spruck type for weak solutions of 

'Av = inC = Q x (0, oo) C 
v = on 8lC = dQ x [0, oo), 
g = v p onQx {0}, ^ ' ' 

v > in C, 

where Q C M™ is a bounded smooth domain, n > 2, and 1 < p < 

For this, we need two nonlinear Liouville theorems for problems in- 
volving the square root of the Laplacian in unbounded domain — one 
in the whole space, another in the half-space. The first one was proved 
by Y.Y. Li, Zhang and Zhu in [18], [17] and Ou in [22J. Its statement is 
the following — and it is equivalent to Theorem II .41 in the Introduction. 

Theorem 6.1. ([IB"]. [XT] . [22]) For n > 2 and 1 < p < 2 tt - 

there exists no weak solution of problem 

t>n+l 



n+1 
n—l ' 




n+1 
n—l '■ 



(6.2) 



We need to prove an analogue nonlinear Liouville type result in- 
volving the square root of —A with Dirichlet boundary value in the 
half-space. This is Theorem 11.51 of the Introduction and Proposition 
16.31 in this section. As we will see, this nonlinear Liouville theorem in 
M." will be first reduced to the one dimensional case M + , by using the 
moving planes method. After this, we prove that there exists no posi- 
tive bounded solution for the nonlinear Neumann boundary problem in 
the quarter which corresponds to the nonlinear Liouville theorem 
involving the square root of —A with Dirichlet boundary value in the 
half-line; see Proposition 16.41 To complete the proof of Theorem 11.51 
we will use the following Liouville theorem in dimension n + 1 = 2. 

Proposition 6.2. ([9]) Suppose that v weakly solves 

-Av>0 mR 2 + , 

^ > on dR 2 + , (6.3) 



v > in Ri 



Then, v is a constant. 



As usual, very strong Liouville theorems (but quite simple to prove) 
hold in low dimensions, but not in higher ones. Compare (16. 3p in 
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low dimensions for supersolutions of the homogeneous linear problem 
with (16. 2p for solutions of a precise nonlinear problem. The proof of 
Proposition 16.21 in [9] compared in an appropriate way the solution v 
with log (| • |). For completeness, we give here an alternative proof. 

Proof of Proposition 16.21 Replacing v by v — inf R 2 v > 0, we may 
assume inf K 2 v = 0. Letting w = 1 — v, we have 

-Aw < in R 2 + , 

$ < on dR 2 + , (6.4) 
w < 1 in R 2 + . 

In addition, sup R 2^ w = 1. Let G C°°(M 2 ) be a function with compact 
support in B 2 r(0), equal to 1 in B R (0), and with |V£r| < ^. Let 

D + R , 2R ■= {(x,y)eR 2 \R< \(x,y)\ <2R,y> 0}. 

Multiplying the first equation in ( 16.41) by w + ^ R , integrating in and 
using the Neumann condition and w + < 1, we see that 

e R \^w + \ 2 <2 / Z R V£ R w + Vw + (6.5) 

1/2 




+ |2 



(6.6) 



This leads, letting R f oo, to L 2 |Vu> + | 2 < oo. As a consequence of 
this, the integral in (16.61) tends to zero as R — > oo. Thus, by ( 16. 5p . 

\Vw + \ 2 = 0. 

Thus, w + is constant, and since sup R ^ w = 1, we conclude w = 1. □ 
Proposition 6.3. Lei n > 2 and 

= = (si, x 2 , ■ • • , x n , y)\x n >0,y> 0}. 
Assume that v is a classical solution of 

(Av = mRlX 1 , 

v = on {x n = 0,y > 0}, 

% = v? on{x n >0,y = 0}, 

w>0 mRl\\ 
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where 1 < p < Then, v depends only on x n and y. 

Proof. We shall follow the steps of [TJ]. Let e n = (0, • ■ • , 0, 1, 0) and 
N — n + 1. Consider the conformal transformation 

z + e n 



z = T{z) 



\z + e r 



and the Kelvin transformation w of v 

w(z) = \z + e n \ N - 2 v{z) = l^l 2 -^). 
Denote B+ /2 (f ) := {z = (x,y) \ \z - \e n \ < \,y > 0}, S+ 2 (f ) : = 

dB+ /2 (f ) n {y > o}, r Q , 1/2 := dB+ 2 (f) n{y = 0}. 

Note that, through T, M"^ 1 = {x n > 0,y > 0} gets mapped into 
the half-ball -B^lif )> the boundary {x n > 0, y = 0} becomes the ball 

ro,i/2, {^n = 0, y > 0} goes to the half-sphere S^df )' anc ^ ^ ne hifinity 
goes to z = 0. 

We see that w satisfies 

(Aw = in (f), 

w = on S+^d), 

5^ = 1^-2)-^^ onrol/2) 

w>0 m5+ 2 (f). 

Since |^|p( n_1 ) _ ( n + 1 ) is nonincreasing in the direction for all i = 
1, - • • , n — 1 (in fact, in any direction orthogonal to the ), the 

moving planes method used as in [9] gives that w is symmetric about all 
the 2 r axis for i — 1, • • • , n — 1. This leads to w = u> (|z'|, z n , y), where 
z! = (zi,-- - , z n -i) an d hence v = v(\x'\,x n ,y). Now, since we may 
perform the Kelvin's transform with respect to any point (— x' Q , —1, 0) 
— and not only with respect to x' = as before — we conclude that 
v = v(x n ,y) as claimed. □ 

Proposition 6.4. Assume that f is a C 1,a function for some a G 
(0,1), such that f > in (0, oo) and f(0) = 0. Let C be a positive 
constant. Then there is no bounded solution of problem 



Av = mR 2 ++ = {x > 0,y > 0}, 

v = on {x = 0, y > 0}, 

% = f{v) on{x>0,y = 0}, 

0<v<C in R 2 ++ . 



(6.8) 



Proof. We use some tools developed in [5]. 

Suppose by contradiction that there is such solution v. First, we 
claim that v (x, 0) — ► as x — ► oo. Suppose by contradiction that there 
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exists a sequence a m — > oo (m — > oo) such that f (a m , 0) — > a > 0. Let 
v m (x,y) := f(x + a m ,y). It is clear that t> m is a solution of (j6.8p in 
f/ m := | x > -a m ,y > 0}. Moreover, f m (0,0) = v(a m ,0) -> a. 

Therefore there exists a subsequence, still denoted by i> m , such that 
v m — > f in C% C (R+) as m — > oo, and f is a solution of 



<9xxf + ^ra^ = ° in K + = y) I y > °} 

< u < C in {y > 0}. 



£ = /(«) > on = 0}, (6.9) 



Notice that 



f(0,0) = a > 0. 

On the other hand, by Proposition 16.21 we know that v is identically 
constant. This is impossible due to the nonlinear Neumann condition, 
since / > in (0, oo) and f(v(0, 0)) = f(a) > 0. We conclude the 
claim, that is, v (x, 0) — > as x — > +oo. 

Note that we can reflect the function v with respect to {x = 0, ?/ > 0}, 
i/) = — v(— x, y) for x < 0, and obtain a bounded harmonic function 
w in all = {y > 0}, since v = on {x = 0, y > 0}. Applying 
interior gradient estimates to the bounded harmonic function v in the 
ball B t (x,t) C R+, we obtain 

\Vv(x,t)\ < G ^°° < j f for all t > i,x > 0. 

On the other hand, by the results of |5J applied to the solution v in 
(or equivalently by the proof of Proposition 13.21 of this paper; note that 
/(0) = 0), we have that |Vu| and \D 2 v \ are bounded in M.\ + n {0 < 
y < 1}. We conclude that |Vf | and |_D 2 t> | are bounded in and 

C 

|Vu(x,t)| < , for alU > 0, x > 0. 

t + 1 

Using interior estimates for harmonic functions as before, but now with 
the partial derivatives of v instead of v, it follows that 

C 

\D 2 v(x,t)\ < , for aUt > 0,x > 0. 

t 2 + 1 

Moreover, we have 

d f\d x v&t)\ 2 -\d y v(x,t)\' 



dx [ 2 
By these facts, we see that the function 



C 

< 



t 3 + l 



THE SQUARE ROOT OF THE LAPLACIAN 37 



is well defined and ^ is also. 

ax 



Using the lim \Vv (x, t)\ = 0, we obtain, for F(v) = f f(s) ds, 



t— *oo 

±m x )+F(v(x,0))} 

+00 

[d xx vd x v - dyvd xy v}(x, t) dt + [f(v)d x v}(x, 0) 

= [d y vd x v + f(v)d x v](x, 0) = 0, 

thanks to the harmonicity of v and the Neumann boundary condition. 
This leads to the Hamiltonian-type identity 

$(•) + F(v(-, 0)) is identically constant in (0, +oo). 

Furthermore, using that lim v(x, 0) = 0, and that lim v(x,y) = 

x— >+oo x— >+oo 

uniformly in compact sets in y (we can prove this by the same previ- 
ous argument leading to lim v(x, 0) = 0), together with the above 

x— >+oo 

bounds for |Vi> (x, y)\ for y large, we deduce 

lim = 0. 

x^+oo 

From all these we obtain 

$(x) + F(v(x, 0)) = 0, for x > 0. 

Since v = and thus d y v = along the y-axis, we see by the definition 
of $(0) that 

i r+oo 

= $(0) + F(v(0, 0)) = $(0) = - / |^| 2 (0, t) dt. 

* Jo 

This implies that d x v = on {x = 0,y > 0}, which contradicts 
Hopf's lemma. Thus, the contradiction means that there is no pos- 
itive bounded solution of the problem. □ 

Before proving Theorems 11.51 and 11.31 let us make some comments. 

Remark 6.5. Theorem 11.51 is still open without the boundedness as- 
sumption on v. 

In this respect, let us give some examples of problems in the quarter 
plane The function v(x,y) = x is an unbounded solution of 

problem 

' -Av = 0, v > in R 2 ++ , 
v = on {x = 0, y > 0}, 

f = on{x>0 )2 / = 0}. 
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This tells us that the result of Proposition 16.21 (which did not require 
boundedness of the solution in the half-plane) does not hold in the 
quarter plane. 

On the other hand, it is clear that v(x,y) = ~ arctan -2^- satisfies 

> for x > 0. Hence, there exists a 



Av 



and —d y v 



|y=0= 



2(l+x 2 ) 



v = 
g > 



bounded harmonic function in the quarter plane such that 

-Av = 0, v > in 

on {x = 0, y > 0}, 
on {x > 0,y = 0}. 

Thus the nonlinear condition |^ = t> p on = 0} is important in 
Theorem 11.51 

Proof of Theorem 11.51 It follows from Propositions 16.31 and 16.41 □ 
Proof of Theorem 11.31 We know by Theorem 15.21 and Proposition 



3Ithat all weak solutions u of (jl.ip belong to C 2 (Q) nC (O). Assume 
by contradiction that the theorem is not true and hence that there is 
a sequence u m of solutions of (11.11) with 

K m = \\u m \\L°°(n) — > oo. 

Since v m = h-ext( Ufri) is £L positive harmonic function in C vanishing 
on 8lC, we have that v m has also K m as maximum in C and that it is 
attained at a point (x m ,0) G x {0}. Let 

and define 

v m (x, y) = K~ 1 v(x m + K m ~ p x, K^ p y), 

We have that ||w m ||L°°(Q m x(o,oo)) < 1 and 

Av m = 0, in C m := Q m x (0, oo) 
v m = on 9fi m x (0, oo), 
on f2 m x {0}, 
in C m . 



x e Vt m ,y > 0. 



v 



dv hi 
> 



(6.10) 



V. 



Notice that 
Let 



v m {0,0) = 1. 



d m = dist(x m , dfl). 
Two cases may occur asm-> oo; either case (a): 

K v ~ x d m — > oo 
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for a subsequence still denoted as before, or case (b): 

Kf^ 1 d m is bounded. 

If case (a) occurs, we have that B K p-i dm (0) = K^ 1 Bd m (0) C Q m and 
that K^ 1 d m — ► oo. By local compactness (Arzela-Ascoli) of bounded 
solutions to (16.101) (recall || , u m ||L oo (n t „) < 1) 5 through a subsequence, we 
obtain a solution v of problem (16.21) in all of = IR n x (0, oo) — note 
that v m (0, 0) = 1 leads to v(0, 0) = 1 and hence v ^ & v > 0. This 
is a contradiction to Theorem 16. 1[ 

Assume now that case (b), K^ x d m is bounded, occurs. Note first 
that since the right-hand side of problem ( 16. ip for v m satisfies \v m \ p = 
v m — -^m; we deduce from the proofs of Proposition 13.11 (iii) and (iv) 
that || Vw m ||xoo(n) < CK'm for a constant C independent of m. Now, 
since u m \ dn = (where u m = v m (-, 0)), we get 

0) < || Vw OT ||i/x>(n)dist(x TO , dVt) < CK^dm. 

We deduce that 

< c < Kv~ x d m 

for some positive constant c. Thus, in this case (b), we may assume 
that, up to a subsequence, 

K v ~ x d m -> a e (0, oo) (6.11) 

for some constant a > 0. 

We deduce that, up to a certain rotation of K n for each index m, since 
we have ( 16. lip . K^ 1 — > oo, d rn — > 0, and that £? K p-i dm (0) is tangent 
to <9f2 m , the domains f2 m converge to the half-space 1R™ = {x n > —a}. 
Thus, through a subsequence of 5 m , we obtain a solution v of problem 
( 16.71) in R"^ 1 = {x n > —a, y > 0} with C bounded by 1 and v > (since 
^m(0, 0) = 1 for all m). This is a contradiction with Theorem 11.51 □ 

Remark 6.6. From Theorem II .31 we have a priori bounds for solutions 
of problem (11.11) with f{u) = u p , 1 < p < ^j-. As a consequence, by 
using blow-up techniques and topological degree theory, one can obtain 
existence of positive solutions for related problems — for instance, for 
nonlinearities f(x, u) of power type, as well as other boundary condi- 
tions. See Gidas-Spruck [Hj for some of these applications when the 
operator is the classical Laplacian. 

7. Symmetry of solutions 

The goal of this section is to prove a symmetry result of Gidas-Ni- 
Nirenberg type for positive solutions of nonlinear problems involving 
the operator A%/2, as stated in Theorem II .61 by using the moving planes 
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method. For this, we work with the equivalent local problem (11.81) and 
derive the following. 

Theorem 7.1. Assume that Q is a bounded smooth domain of R" 
which is convex in the x\ direction and symmetric with respect to the 
hyperplane {xi = 0}. let f be Lipschitz continuous and let v G C 2 (C) be 
a solution of (11.81) . where C = Q x (0, +oo). Then, v is symmetric with 
respect to X\, i.e., v(—Xi,x',y) = v(xi,x',y) for all (—Xi,x',y) G C. 
In addition, J^- < for X\ > 0. 

Proof of Theorems 11.61 and 17. 1L It suffices to prove Theorem 17.11 
From it, Theorem 11.61 follows immediately. 

Let ) G f2 and A > 0. Consider the sets 

Sa = {(x\,x) G Q | x\ > A} and T\ = {(xi,x) G ft \ X\ — A}. 

For x G Sa, define X\ = (2A — Xi,x'). By the hypotheses on the 
domain Q we see that 

{x\ | x G S A } C Vt. 

Recall that v G C 2 (C) is a solution of 

Af = in C = fix(0,oo), 

v = on d L C = dQ x [0, oo), 

£ = /(«) on fix{0}, 

■u > in C. 

For G Sa x [0, oo), let us define 

v x(x, y) = v(x x , y) = v(2X - x l} x', y) 

and 

w\(x,y) = Oa - v)(x,y). 

Note that v \ satisfies 

Av\ = in S A x (0, oo), 
v x >0 on (dttf) S A ) x (0,oo), 

!& = fM on S A x{0}. 

Thus, since <9£a = (dQ fl £a) U Ta and w\ = on T\, we have that w\ 
satisfies 

Au'a = in S A x (0, oo), 

^a>0 on (<9£ A ) x (0,oo), (7.1) 



where 



^ + c a (x)«;a = on E A x{0}, 
/Oa) - /(«) 



c A (x, 0) 



Note that c A (x,0) G L°°(S A ). 
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Let A* = sup{A | Sa 7^ 0} and let e > be a small number. If 
A G (A* — e, A*), then Ea has small measure and we have, by part (ii) 
of Proposition 14.41 (applied with Q replied with E A ), that 

w\>0 in E A x (0, oo). 

Note here that Ea is not a smooth domain but that part (ii) of Propo- 
sition 14.41 does not require smoothness of the domain. By the strong 
maximum principle, Lemma 14.21 for problem (17. ip we see that W\ is 
identically equal to zero or strictly positive in Ea x (0, oo). Since A > 0, 
we have w\ > in (dQ n <9Ea) x (0, oo), and hence we conclude that 
w\ > in E A x (0, oo). 

Let A = inf{A > | w\ > in E A x (0, oo)}. We are going to 
prove that A = 0. Suppose that A > by contradiction. First, by 
continuity, we have W\ > in Ea x (0, oo). Then, as before, we deduce 
w\ > in Ea x (0, oo). Next, let 5 > be a constant and K C Ea 
be a compact set such that |Ea \ K\ < 5/2. We have W\ (-, 0) > r] > 
in K for some constant rj, since K is compact. Thus, we obtain that 
w\ - £ (-, 0) > in K and that |Sa _ £ \K\ < 5 for e small enough. 

Now we apply again part (ii) of Proposition 14.41 in Ea _ £ x (0, oo) 
to the function w\ - £ - We know that u>a o _ £ (-,0) > in K, and hence 
{^Ao-e < 0} C Ea _ £ \ K, which has measure at most 5. We take 5 to 
be the constant of part (ii) of Proposition 14.41 We deduce that 

wxo-e > in Sa - £ x (0, oo). 

This is a contradiction to the definition of Ao- Thus, Ao = 0. 
We have proved, letting A j A = that 

v(—xi, x', y) > v(xi, x', y) in (Q fl {x\ > 0}) x (0, oo) 

and, since w\ = on T\, 

cLf = — ~7T^- < for x\ > 0, 
2dx 1 

by Hopf's lemma. Finally replacing x\ by —x±, we deduce the desired 
symmetry v(— xi, x', y) = v(xi,x', y). □ 
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